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ABSTRACT 

The polarised spectrum of stellar radiation encodes valuable information on the condi- 
tions of stellar atmospheres and the magnetic fields that permeate them. In this paper, 
we give explicit expressions to estimate the magnetic field vector and its associated 
error from the observed Stokes parameters. We study the solar case where specific 
intensities are observed and then the stellar case, where we receive the polarised flux. 
In this second case, we concentrate on the explicit expression for the case of a slow 
rotator with a dipolar magnetic field geometry. Moreover, we also give explicit formu- 
lae to retrieve the magnetic field vector from the LSD profiles without assuming mean 
values for the LSD artificial spectral line. The formulae have been obtained assuming 
that the spectral lines can be described in the weak field regime and using a maximum 
likelihood approach. The errors are recovered by means of the hermitian matrix. The 
bias of the estimators are analysed in depth. 
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1 INTRODUCTION 

The most precise measurements of stellar magnetic fields 
. are based on the observation and the interpretation of 
polarisation in spectral lines. Most of the magnetic fields 
detected in stars at different phases of evolution have 
been found by means of the Zeeman effect, which gener- 
ates linear and circ ular polarisation i n the presence of a 



magnetic field (e.g., Wade et al. 2000; Bagnulo et al. 2002; 



Jordan. Werner fe O'Toolel 120051: Aznar Cuadrado et al 



20041 : lO'Toole et al.ll2005l : ISilvester et alj|2009l : iLeone et a! 



20111 ). Here we present analytical expressions for infer- 



ring the magnetic field vector from the observed Stokes 
profiles induced by the Zeeman effect in the weak field 
approximation. 

The weak field approximation is broadly applied for the 
inference of solar and stellar magnetic fields from the obser- 
vation of the Stokes profiles. It is an analytical solution to 
the radiative transfer equation whose main basic assump- 
tion is that the magnetic field is sufficiently weak throughout 
the whole region of the atmosphere where a spectral line is 
formed. Although simple, this approximation is very useful 
since non-magnetic mechanisms usually dominate the shape 
of spectral lines, in both the solar and the stellar case. 

For example, in the quiet Sun, which occupies the 
vast majority of the solar surface (far from the sunspots 
that harbour very strong fields), the spectral lines are well 
described by the weak field approximation. This has allowed 
the success of many synoptic magnetographs like those 



of Big Bear l|Spirock et all l200ll : IVarsikl 1 19951) . It is even 
used to produce modern vector mag netograms like those 
obtai ned with the IMaX instrument (iMartmez Pillet et al.l 
l201ll ) onboard the Sunrise balloon i|Solanki et al.l |201Ch . 
The weak field approximation has also been used 
to diagnose the chromosphere dMerenda et all 120061 : 
Asensio Ramos. Truiillo Bueno fe Landi Degl'Innocentil 
2008) due to the enhanced thermal width of the spectral 
lines. 

In night-time spectro-polarimetry, the weak field ap- 
proximation is at the base of t he least-squares deconvolu- 
tion (LSD; iDonati et all 1 19971 ). the most successful tech- 
nique used to detect and measure magnetic fields in solar- 
type stars or other stars in which the polarisation signal 
per spectral line is well below the noise level. Other works 
have used this approximation to diagnose magnetic fields 
in a large variety of stellar objects. A limited selection in- 
cludes ^s£mej^cejrt_woj^cs^^ 

ulae dJordan. Werner fe O'Toolel 120051: ILeone et all feoill ). 
white dwarfs l|Aznar Cuadrado et al.ll2004 | ). pulsating star s 
|Silvester et al.ll2009l), hot subdwa rfs (jO'Toole et alj|2005h . 
Ap an d Bp star s (Wade et al. 2000) and chemically peculiar 
stars (Bag nulo et al.ll2002l ). 



2 THE WEAK FIELD APPROXIMATION 

The weak field approximation is an analytical solution to the 
radiative transfer equation. The fundamental assumption is 
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that the magnetic field vector is constant and its intensity 
is sufficiently weak throughout the whole region of the at- 
mosphere where an spectral line is formed. Additionally, the 
line-of-sight velocity and any broadening mechanism have 
to be constant with height in the line formation region. As 
a consequence, the magnetic field can be considered as a 
perturbation to the zero-field case. Quantitatively, the ap- 
proximation holds whenever A As <C A Ad, where A As is 
the Zeeman splitting and AA_d is the dominant broadening 
mechanism (thermal, rotation, etc.). From its definition, it 
is clear that the weak field regime occurs at different field 
strengths for different spectral lines (depending on the sen- 
sitivity to the magnetic field, the local temperature and the 
atomic mass) and for different stellar objects (depending 
on any non-magnetic broadening mechanism for the spec- 
tral lines, the rotation being the most efficient mechanism 
in cool stars). 

To first order in A As, the intensity profile of a spectral 
line formed in a weak magnetic field is insensitive to the 
magnetic field. In other words, it fulfils the transfer equation 
in the absence of a magnetic field. At this first order, the 
circular polarisation profile, i. e., the Stokes V profile, for a 
given spectral line has the following expression: 

dl(x) 



V{x) = -CAgB\\ 



Ox 



(1) 



The symbol Bu stands for the longitudinal component of 
the magnetic field, i. e., Bu = Bcos9b, where 9b is the 
inclination of the magnetic field with respect to the ob- 
server's line of sight (f2; see Fig. [TJ and B is the mag- 
netic field intensity. The symbol g represents the effective 
Lande factor of the line, which quantifies the magnetic 
sensitivity of the line. The effective Lande factor only de- 
pends on the quantum numbers of the transition (see, e.g., 
lLandi Degl'Innocenti fe Landolfil [2004). Equation [1] is ex- 
pressed in terms of the generic wavelength variable x. If 
x represents the wavelength A, then the parameter A equals 
the central wavelength of the line, Aq. However, it is cus- 
tomary in stellar spectropolarimetry to consider x as the 
velocity in Doppler units. In such a case, we have A = cAo, 
with c the speed of light. The constant C = 4.67 x 10~ 13 
G- 1 A" 1 . 

At first order in A As, both Stokes Q and U are zero. In 
order to obtain an expression for the Stokes profiles charac- 
terizing linear polarisation we have to expand the radiative 
transfer equation to second order in AAs and assume that 
the spectral line is not saturated. Under these assumptions, 
the general formulae for the Stokes Q and U are: 



Q(x) = ~^-A 2 GB 2 ± cos20s^^, 
4 ~ ox 2 

U{x) = „<±. A 2 GB 2 ± sm2d? B d2l{x) 



dx 2 



(2) 



The symbol G plays the role of the effective Lande factor 
for linear polarisation and quantifies the sensitivity of linear 
polarisation to the magnetic field. Again, it is only a func- 
tion of the quantum numbers of the transition (see, e.g., 
lLandi Degl'Innocenti fc Landolfil [20041 ) . The symbol B± = 
B sin 9b is the component of the magnetic field perpendicu- 
lar to the line of sight. The angle 4>b represents the azimuth 
of the magnetic field vector with respect to an arbitrary ref- 
erence direction (e a , ej in Fig.[T]show the coordinates chosen 
for the reference of Q > 0). 



When observing the polarised light in resolved sources 
like the Sun, we detect specific intensities. However, when 
observing unresolved stars, the incoming polarised radiation 
is obtained as an integration in the plane of the sky of the 
individual Stokes parameters at each point of the stellar sur- 
face. As a consequence, the specific values of the Stokes flux 
vector depends on the surface distribution of the magnetic 
field, on the centre-to-limb variation (CLV) of the radiation, 
and the Doppler effect due to stellar rotation. In this pa- 
per we neglect rotation, the derived expressions being valid 
for those cases where the broadening mechanisms dominate 
over rotation. In the Sun, since we observe local profiles, the 
main broadening mechanism is typically thermal. In gen- 
eral, the polarimetric signals of 99% of the solar surface, the 
so-called quiet Sun, can be explained in terms of the weak 
field regime in most spectral lines in the visible and the 
near-IR. For non-resolved objects, the applicability depends 
strongly on the actual rotational velocity as compared with 
the temperature, the observed spectral line, the spectral res- 
olution, the magnetic field organisation, etc. The thermal 
broadening depends on the square root of the ratio between 
the temperature and the atomic weight of the atom. There- 
fore, in general, the larger the temperature and the lightest 
the atom, the larger the allowed rotational velocity. For low 
resolution spectrographs (R ~ 2000), basically the only re- 
maining spectral lines (that are not wiped out by the lack 
of spectral resolution) are from hydrogen. For these lines, 
for a temperature of T ~ 10000 K, the maximum line of 
sight velocity is 10 — 20 km s _1 . However, when observing 
metal lines at spectral resolutions as high as R ~ 60000, the 
allowed velocities are one order of magnitude lower. 

It is customary to introduce a parametrised for m of the 
CLV . In this wo rk, we assume a quadratic form (e.g.. IClaretl 
l200d : ICoxll2000r i, which gives a good balance between the 
quality of the CLV and the simplicity of the analytical ex- 
pressions presented in this paper for the polarised fluxes: 

I(x,fi) = Io(x) (l — u — v + Uii + vp 2 ) = I (x)f(n), (3) 

where Io(x) is the intensity profile at disc centre (p = 1). 
The parameters u and v have values between and 1 and 
are supposed to be constant along the spectral line (although 
they can vary from line to line). Note that the values of u 
and v have to fulfil the condition that I(x) > 0. The CLV 
is given in terms of p = cos O, where Q is the astrocentric 
angle between the normal to a point in the stellar surface and 
the line of sight. Using this law for the CLV (and neglecting 
rotation) the flux for Stokes / is: 



Fi{x) = Ioix) / dE'fifi) 



(4) 



where the integral is computed over the visible surface in 
the plane of the sky E'. For an arbitrary function h(p, a) of 
the polar coordinates of the stellar surface p and a, we have 



// 



r 1 p2tt 

dT,'h(p,a) = / pdp / dah{p 7 a). 
Jo Jo 



(5) 



Note that the variable p = sin <d, which means that p = 
\J\ — p 2 . Plugging the quadratic expression considered for 
the CLV, the final closed expression for the integrated Stokes 
/ is found to be: 



Fi(x) = iv Ioix) 



6 — 2u — 3v 



(6) 
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Figure 1. Geometry of the stellar model. The symbol f2 represents the line of sight. The vector B displays the magnetic field vector in 
the stellar surface. Its geometry is described in terms of the inclination with respect to the line of sight 9 b and the azimuthal angle <j>g 
(the axis e a being the reference for the azimuth). The position of a point P in the stellar surface (£) is defined only by the astroccntric 
angle 0. Its projection (P') on the plane of the sky (S') is represented in cylindrical coordinates by the modulus p and the angle a, 
which is referring to the e a and e*(, axis. 



Integrating Eqs. JTJ and ([2]) over the visible surface, we end 
up with the following expressions for the polarised flux: 



Fv{x) = -GAg{Bu) 



dx 



-A G{B± cos20 fl )- 



9a; 2 



C 2 .2^ /D 2 _._ OJL ,d 2 Fi(x) 



F v {x) = —-—A G{Bx sui20b) 



with the definitions: 



<5||) 



{B± cos 20b) = 
(B 2 ± sin 20 B ) = 



6 


tt(6- 


2u 


-3v) 




6 




tt(6- 


2u 


-3v) 




6 




tt(6- 


2u 


-3v) 



J dS' B l{ f(fi) 

dE' B\ cos20b/(ai) 
dS' B\ sin 2(j> B /(/i). 



(8) 



Therefore, the weak-field approximation for the integrated 
polarised flux remains formally the same but the compo- 
nents of the magnetic field Bu and B± appear weighted 
by the CLV law. In general, the components of the mag- 
netic field may depend on the position on the stellar sur- 
face in a complicated manner, so that the previous inte- 
grals might not have closed expressions. In the simple case 
that the magnetic field is constant across the stellar surface, 
we recover {Bu) = By, (-B 2 _cos20b) = B\cos2<f)B, and 
{Bj_ sin 20b) = B\ sin 20b- One of the simplest non-trivial 
configuration we can consider explicitly is that of a dipolar 
field. The magnetic field vector at each surface point is then: 



B(r) 



Hd 
2 



e — 3(e • r)r] 



(9) 



where the unitary vector e defines the orientation of the 
dipole and the unit vector r indicates positions on the stel- 
lar surface. The quantity Hd represents the magnetic field 
strength at the poles of the dipolar field. From Eq. [7] and 
Eq.[9]and after some algebra, it is possible to write the flux 



of the Stokes parameters as: 

, . 1 15 + u dTi(x) 

Fv[x) = --CA 9 6 _ 2u _ 3v H,^ r 



Tu(x) = 



C 2 ,2^420 - 68w - 105w rr2 nJ _ d 2 Ti{x) 



4480 
'4480" 



-A G- 



6 — 2u — 3v 



-H± cos 20d- 



dx 2 



A 2 G^- 



lu-W5v u2 . d 2 Ti{x) 
-H ± sin 20 d - 



2u — 3v 



dx 2 



(10) 



where Hu = HdCosOd and H± = H d sin 8d, Qd being the 
inclination of the axis of the dipole whi le <j>d is the azimuth o f 
the dipole (for a similar derivation see lLandolfi et aHll993T ). 

The previous formalism has demonstrated that the 
weak-field approximation leads to formally the same ex- 
pressions either if we consider specific intensities (that is 
applied whenever spatial resolution is available) or if we 
consider fluxes (whenever the object of interest cannot be 
resolved). In other words, the observed circular (linear) po- 
larised spectrum is proportional to the first (second) deriva- 
tive of the observed intensity through the longitudinal (or- 
thogonal) component of the magnetic field. The only differ- 
ence resides on the exact definition of the components of the 
field which, in the spatially unresolved case, are an average 
over the stellar surface weighted by the CLV. Therefore, for 
the sake of simplicity, it is possible to combine all expressions 
on the following general ones: 



V = -CB\\T' 

Q = -C 2 B\ cos 201" 



U 



-C 2 B\ sin 202" 



(11) 



The meaning of the newly defined variables is summarized 
in Tab. Q] 



3 ESTIMATION OF THE MAGNETIC FIELD 
VECTOR 

Once the model is set, our aim is to infer the magnetic 
field vector parametrised in terms of Bu, B±, and from 
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Table 1. Variables in Eq.lTTI 



Variable 



Resolved 



Dipole 



I 
V 
Q 
u 

X 1 
X" 

By 



i 

v 
Q 

OX 

B 



B 



1 15+u 



A dZiM. 

10 6-2u-3v y dx 
420-68u-105« ^2 ( ^ d 2 .7 r J<>) 



-3d 

Hx 

4>d 

d 



observations of a set of spectral lines (X 1 , Q* ,lf ,V*) with 
i = l... niines. Assuming that the weak field approximation 
can be applied for the set of observed spectral lines and 
that observations are corrupted with uncorrelated Gaussian 
noise, we can use a least-squares estimator (maximum likeli- 
hood) to retrieve the magnetic field vector (see the Appendix 
for more details). The \ 2 merit function is defined as: 



(vj - v;'' mod ) 2 



E 



(i ji 7 ji;mod\2 

2— 1 (rr* "I 2 ' 



(12) 



where the subindex j indicates the position along the coordi- 
nate x. The label mod refer to the model for Stokes profiles, 
and the superindex i labels the spectral lines. The previ- 
ous merit functions consider the quite general case in which 
the standard deviation of the corrupting Gaussian noise is 
different for Stokes Q, U, and V and for each wavelength 
point Xj. However, in practice, we have the simpler situa- 
tion in which a 



Vj 



j. Moreover, although 
the number of photons arriving in the line cores are smaller 
than in the far wings, we make the approximation that the 
noise variance is wavelength-independent. Furthermore, we 
also consider that it is independent of the considered line, 
so that cr| = a. These simplifications lead to less cluttered 
expressions for the inferred parameters. However, we point 
out that the general expressions that emerge from the opti- 
mization of Eq. (|12p can be found in the Appendix. 

In order to infer a certain parameter, we have to find 
the global minimum of the \ 2 ■ This is trivially obtained 
by solving the non-linear system of equations obtained by 
setting the derivatives of the \ 2 function with respect to 
that parameter to zero. By so doing (see Appendix), we can 
obtain the expression of the magnetic field vector in terms 
of the observables: 



Bn = 



i Y.,,v r. 
CT.^')) 2 



(£« C^"5) a + (E«^"j: 



c* 2 



1 Znup 1 , 

- arctan — — - — : + mo , 

2 E, 3 e;z"3 



(13) 



and the derived quantities: 
Bx 



tan 6 — 



B = ^B\ + B- 



(14) 



The phase shift (j>o is used to set the correct quadrant for 
the azimuth and depends on the sign of the numerator 



j LtjZ") an d the denominator £2 



lows. If £2 / then: 



Ei, QV") as fol- 



In the case £5 




if it 5s and £2 > 
if U < and £2 > 
if £2 < 



\ if U > 



3tt 
4 



if it < 



(15) 



(16) 



If both £2 and it are zero, the angle <f> is obviously undefined. 

The estimate d errors can be com puted from the covari- 
ance matrix (e.g.. [Press et all 1 19861 ) . In the simple case we 
consider of equal wavelength-independent standard devia- 
tion for all spectral lines, the covariance matrix is diagonal 
(no correlation between the parameters), and the errors at 
a confidence level of 68.3% (one sigma) are expressed as: 

5B U 1 
SBx 



a 



2C 2 BxJJ2 i] (l" i j ) 2 



2C7 2 BWE«(2'1) 2 



= ± 



B\SB 2 +B\8B\ 
B 2 +B 2 



SB = ± 



\ 



B 2 SB 2 +B 2 ± 5B 2 
Bl +B 2 



(17) 



4 BIAS OF THE MAXIMUM LIKELIHOOD 
ESTIMATOR 

It is widely known that all maximum likelihood estimators 
may suffer from biases. The bias is the difference between the 
value of the estimator and the true value of the parameter. 
Each individual parameter has to be studied separately us- 
ing analytical/numerical simulations to understand to which 
extent the estimation of B\\ , Bx, and 4> are subjected to bias. 
For simplicity the simulations we present in the following re- 
fer to the resolved case in which we use specific intensities. 
However, the behaviour of the estimator is the same also for 
the stellar case. 

In order to carry out the simulations, we focus on the Fe 
I line with central wavelength Ao = 5250.2 A. This spectral 
line is produced by the transition 5 Do — 7 D\. The values of 
the effective Lande factor for circular polarisation is g — 3, 
while it is G = 9 for linear polarisation. We consider a con- 
stant magnetic field strength of B — 500 G, which is suffi- 
ciently weak so that the weak-field approximation can still 
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Figure 2. Bias of the maximum likelihood estimator of the magnetic field in the weak field approximation. Different color lines represent 
different values of the noise level added to the synthetic profiles of the Fel line at 5250.2 A . The color code is the following: yellow 
represents the inversion of the profiles without noise. Orange displays the results for a noise level of f0 — 5 f c , pink for 5 X f0 — 5 I c , red for 
10~ 4 l c , brown for 5 X 10~ 4 I c , dark green for 10 -3 f c , light green for 1.5 X 10~ 3 f c , light blue for 2 X 10~ 3 l c , dark blue for 5 X 10~ 3 
f c , light violet for 10 — 2 I C) and dark violet for 5 X 10 — 2 I c . The l c is the continuum intensity. 



be considered. The inclination and the azimuth of the field 
are set to vary uniformly between and 180°. Letting the 
azimuth vary in this interval, we avoid the 180° ambiguity 
in the azimuth present in the radiative transfer equation. 
For simplification, we consider a Gaussian intensity profile 
in the form: 



7(A) 



d c exp 



(A-Ao) 2 



2w 2 



(18) 



where d c = 1/2 and w = 0.05 A. The parameters for the 
Gaussian profile have been fixed to fit a solar observation 
obtained with the IMaX instrument. From the IMaX obser- 
vational capabilities, it can be verified that the 5250.2 A line 
is in the weak-field regime up to ~ 1 kG (with the line as- 
sumed to be in local thermodynamical equilibrium in a quiet 
Sun model atmosphere). We have synthesized 500 profiles for 
different combinations of the inclination and azimuth and 
different noise realisations. The noise added has a Gaussian 
statistics and we consider the effect of different standard de- 
viations. We use Eqs. (|13p to compute the inferred values of 
the parameters. Since we repeat the experiment for differ- 
ent realisations of the noise, we end up with a distribution 
of values for each parameter. We adopt the median value as 
the estimation of the parameter (the percentile 50 P50, i. e. 
, the value of the parameter that contains 50% of the area 



of the distribution). To quantify the dispersion produced by 
noise we use the percentiles 16 and 84 (which encompass one 
standard deviation around the estimated value). 

Figure [2] displays the inversion of the synthetic profiles 
with different values of the standard deviation of the noise, 
from to 5 x 10~ 2 in units of the continuum intensity, I c . 
It is clear from the upper and lower left panels that both 
the longitudinal component of the magnetic field and the 
azimuthal angle are unbiased quantities. This means that 
the estimated value statistically coincides with the original 
one. As expected, the dispersion grows with the increasing 
noise level. 

Contrary to By and (j>B (note that the notation for 
the magnetic field vector is the one associated to resolved 
sources), the transversal component of the magnetic field 
and hence the inclination angle presents a non-zero bias. 
Except for the case in which there is no added noise, small 
transversal components of the field (smaller or at the level 
of the noise amplitude) are overestimated. The fundamental 
reason is that the expression for B± in Eq. (|13[) is not ro- 
bust against noise. It can be verified that if and U] are at 
the noise level (can be described as Gaussian distributions 
with zero mean and variance a 2 ), B± is a random variable 
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following the probability distribution: 



using: 



P(B X ) = ^4 L B 3 ± e X p 



2<r 2 



(19) 



The value of this distribution for a given percentile c fulfils: 

!/4 _ 



-21n(l - c) 

£ y (Z"5) 2 



(20) 



The percentile 50 (c = 0.5) correctly captures the value of 
the bias at small values of B±. Once this value is computed, 
if the inferred value of B±_ is similar to this value, one should 
be aware that the correct value of B± might be smaller. If 
this is not taken into account, the estimated inclination of 
the field is larger than the real one and artificially horizontal 
fields might be inferred. Note also that in the stellar case, 
the inferred inclination of the dipole axis would be more 
inclined than the real one. 



5 THE PARTICULAR CASE OF THE 
LEAST-SQUARES-DECONVOLUTION 
PROFILE 

Most detections of faint signals in stellar atmospheres have 
been evidenced by adding many spectral lines. The polari- 
metric signal per spectral resolution clement is known to be 
well below the noise level, so that line addition is a must to 
fight against photon noise. The most widely used and suc- 
cessful technique that combines the information of many 
spectral li nes is the Least-S quares Deconvolution (LSD) 
technique (|Donati et al.| [l997). The equations presented in 
this paper allow us to retrieve the magnetic field vector from 
the polarised stellar spectra taking into account many spec- 
tral lines. Thanks to it, we can rewrite the equations to 
directly extract information from the LSD profiles. 

The LSD technique is fundamentally based on the appli- 
cation of the weak-field approximation and on the assump- 
tion of a constant CLV for all spectral lines (in our case, u 
and v are set constant). This means that all Stokes profiles 
for each spectral line can be computed with a single spectral 
profile whose proportionalit y constant changes from line to 
line (see lDonati et al.lll997h : 



T = rfl 
V = T) i A i g i V 
Q l =rf(k i fG i Q 
It = rf{K i f(TU, 



(21) 
(22) 
(23) 
(24) 



where rf is the line depth and X, Q, U, and V are the LSD 
profiles that can be c omputed usin g a le ast-squares pro- 
cedure as explained in [Ponati et al.l (1 19971) . We refer to 
iKochukhov. M akaganiuk & Piskunov (201^) for an 
in-depth analysis of the assumptions and potential 
problems of LSD. From the previous equations, it is easy 
to show that it is possible to infer the magnetic field vector 



B ± = 



1 J2j Vj 



CK 



9 2 Z 3 



C 2 K' 



1 J2iUj-g^T , , 
- arctan _ 2± + <po, 



(25) 



where K = 1 and K 1 — 1/4 for the case of an unresolved 
star with constant magnetic field and 



K = 



15 + u 



10 6 - 2u - 3v ' 



K' = 



1 420 - 68u - 105w 



4480 



6 — 2u — 3v 



(26) 

for the stellar dipole. Note that it is not necessary to assume 
averaged atomic parameters for the LSD profile, treating it 
as a mean spectral line. In our case, the estimated magnetic 
field vector only depends on the observables, the atomic pa- 
rameters of each observed spectral line (which is necessary 
to compute the LSD profile), and the assumed CLV coeffi- 
cients. 



6 ILLUSTRATIVE EXAMPLES 

The inference power of the expressions developed in 
the previous sections are illustrated with the aid of 
two different examples. The first one consists of a 
simulated stellar dipole using a M ilne-Eddington (e.g., 
Landi Dcgl'Innoccnti & Landolfi 2004) atmosphere and the 
second one is a particularly interesting observational exam- 
ple in which we can illustrate the effect of the bias in the 
transverse component of the magnetic field with spatial res- 
olution. 

Considering the stellar dipole, we simulate a Milne- 
Eddington atmosphere with a source function that varies 
as: 



S(t)=So(1 + Pt). 



(27) 



For simplicity, we choose /3 = 1 and consider a static atmo- 
sphere at all points of the stellar surface. We assume a spec- 
tral line centred at A = 5000 A with a Doppler broadening of 
0.04 A . Both the circular and linear effective Lande factors 
are equal to 1. The CLV variation of the Milne-Eddington 
atmosphere (that we have forced to be wavelength indepen- 
dent) gives u = 2/3. We integrate the Stokes signals on the 
visible stellar surface and we add Gaussian noise to the final 
integrated flux with a standard deviation of 5 x 10 -5 in units 
of the intensity flux at the continuum, J- c . 

Figure [3] shows the flux of the simulated Stokes param- 
eters coming from a dipole field with the following parame- 
ters: H d = 1500 G, d = 80°, and <j> d = 25°. The black lines 
represent the synthetic fluxes, and the rhombs the synthetic 
observations with added noise. We invert the noisy Stokes 
fluxes using Eqs.[H[T4]and[T7]and obtain H\\ = 265.5 ± 2.4 
G, H± = 1501.6 ± 25.7 G, </> d = 27.2° ± 2.0, and the derived 
quantities H d = 1524.9 ± 25.3 G and d = 80.0° ± 0.2. All 
quantities are nicely recovered. In fact, the bias estimations 
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Figure 3. Inversion of a synthetic stellar dipole. The black lines display the synthetic flux of a dipole with a strength H d =1500 G, an 
inclination 8 d = 80° with respect to the line of sight and an azimth of 9 d = 25°. The rhombs are the synthetic polarised flux including 
an additive Gaussian noise with a standard deviation of 5 X 10 — 5 T c - The blue lines are the fluxes inferred from the inversion using Eq. 
13. 



for the perpendicular component H± for a the percentiles 
16, 50, and 84 (i. e. containing one sigma probability) are 
302.0, 426.5, and 543.8 G, respectively. This means that the 
computed H± is reliable. Now, we assume a weak dipole 
with H d = 100 G, 9 d = 20°, and 4> d = 25°, and the same 
noise level. The inferred parameters are H\\ = 99.7 ± 2.3 G, 
H± = 333.9T 109.0 G, 4> d = 39.8° ±37.4, H d = 348.3T 104.5 
G, and 9 d = 73.4° ± 5.1. In principle, both the longitudinal 
and the transverse components (as well as the inclination) 
should be well recovered but the azimuth remains undeter- 
mined. However, the bias of H± for the percentiles 16, 50, 
and 84 are 293.2, 414.1, and 528.0 G, respectively. There- 
fore, the inferred perpendicular component of the dipole is 
consistent with a bias. This implies that the inferred value, 
although having a small (Gaussian) error, it has to be con- 
sidered an upper limit. In this case, we know that the per- 
pendicular component of the dipole is very small (34 G). 
Consequently, we overestimate both the perpendicular com- 
ponent and the intrinsic strength of the dipole. Additionally, 
the dipole appears to be much more inclined that in reality. 

The second experiment consists of filter-polarimetric 
data observed with the IMaX instrument on the Sunrise 
mission. This polarimeter observes the 5250.2 Fel line (for 
which we have carried out the bias experiment in Section 
4). The data set consists of the four Stokes parameters ob- 
served at the quietest areas of the solar disc centre at a spa- 
tial cutoff frequency of about 0.15-0.18" (~ 120 km in the 
solar surface, the best one at the moment for instruments 
with polarimetric capabilities). The noise level in circular 



and linear polarization is 10 -3 in units of the continuum 
intensity, I c . The left panel of Fig. [4] shows the continuum 
intensity image, with brighter areas associated to granular 
regions, where the plasma ascends to the photosphere. Dark 
areas are the intergranular lanes, where the motions are pref- 
erentially downflowing. The right panel of this same figure 
displays the estimated bias for the transversal component of 
the magnetic field using Eq. [20] for the percentile c = 0.5. 
The bias has a very particular spatial distribution which 
mimics reversed granulation (bright areas become dark and 
viceversa). This particular spectral line is strongly sensitive 
to the temperature, becoming very deep and narrow in in- 
tergranular lanes and less deep and broad in granules. The 
dependence of Eq.[20]on the second derivative of the inten- 
sity profile, which gives an idea of the width of the spectral 
line, is the reason why the bias is larger in intergranules and 
less important in granules. 

Left panel of Fig. displays the inferred longitudi- 
nal component of the magnetic field. As can be seen, the 
noisy background has values around zero, consistent with 
the fact that the estimator of this quantity is unbiased. 
The central panel shows the inferred transversal compo- 
nent of the magnetic field. The noisy background is now 
filled by magnetic fields that have quite intense values, il- 
lustrating the network pattern of the bias. If this bias for 
the longitudinal field is not appropriately account for, this 
might lead to an artificial excess of inclined magnetic fields. 
This effect could be affecting some of the recent magnetic 
field inferences in the quiet Sun observed with Hinode (see, 
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Figure 5. Left (centre): inferred value of the longitudinal (transverse) magnetic field for the 5250.2 A line. Right: transverse field 
corrected from the bias due to pure noise. Central and right panels share the same colorbar. 



e.g JOrozco Suarez et al. l l2007l ; iLites et al.ll2008l ; ISheminoval 
bOOallshikawa & Tsunota 2009, 20ld. l201ll)~ 



Since we can characterize the bias by its median value, 
it is possible to distinguish real signals from false signals pro- 
duced by the presence of noise. The way we proceed is as 
follows. We compute a conservative upper limit for a trust- 
ful field as the bias for the percentile 84 (c = 0.84 in Eq. 
I20[) . This value changes from pixel to pixel. Then, we force 
B± — in those places where the inferred B± is smaller than 
the bias. The result is represented in the right panel of Fig. [5] 
Now, the real signals (coming from pixels with linear polar- 
ization clearly above the noise level) are much more evident 
and most of the background has disappeared. Note that we 
have only removed those signals that were produced by the 
presence of noise (the expected value of B± = should be 
zero). However, if B± ^ 0, the bias can still be important if 
the noise level is high, as shown in Fig. [2] 



papers is that we give explicit formulae for the components 
of the magnetic field vector in terms of the observables. The 
formulae are general and hold for specific intensities, for in- 
tegrated fluxes and are slightly modified for LSD profiles. 
In the particular case of a stellar dipole, the orientation of 
the dipole axis and its intensity can be recovered from one 
observation of the full Stokes vector. 

We have also studied the bias of this maximum like- 
lihood estimator. The longitudinal magnetic field and the 
azimuthal angle are unbiased quantities. However, the 
transversal component of the magnetic field and hence the 
inclination of the field are overestimated in the presence of 
noise. We derive the estimated value of the perpendicular 
component of the magnetic field in the case that there is 
no linear polarization signal above the noise (the bias for 
B±=0). We propose to evaluate this bias prior to the infer- 
ence of the perpendicular component of the field. One should 
be very cautious when the inferred B± is of the order of the 
bias. 



7 CONCLUSIONS 

We have shown that the weak field approximation (the 
Stokes parameters are proportional to the first and second 
order derivative of the intensity) also holds for observed stel- 
lar fluxes in the case of slow rotators. We have used a max- 
imum likelihood estimator to infer the magnetic field vector 
from the observed Stokes profiles. The main result of this 
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APPENDIX A: DERIVATION OF THE MAXIMUM LIKELIHOOD SOLUTION OF THE STELLAR 
MAGNETIC FIELD 

We start from the general weak-field equations that hold both for the resolved (solar) case and for the integrated stellar dipole 
(Eq. 11): 

V = -CB\\T' 

Q = -C 2 B\ cos 201" 



U 



-C 2 B\ sin 201" 



(Al) 



where the explicit expressions for V, Q,U,X', and X" , for each case can be found in Table [T] in the main text. We denote 
the Stokes vector by S = [I, Q,U, V]. Assuming that the difference between the data and the model follows a Gaussian 
distribution, the likelihood (the probability distribution of the data given the parameters) can we written as: 



£ = 



nnnM fc ) 2 p /2 exp 



i=l j=l fc=2 



J ik 



i\2 



2(^) 2 



(A2) 



where the index i refers to the spectral line, j refers to the wavelength, and k to the element of the Stokes vector. The 
abbreviation "mod" sta 
the log-likelihood, ln£: 



abbreviation "mod" stands for the model. The symbol a 2 stands for the variance of S - 5 mod . Taking the logarithm, we build 



3nmx 



i=l j=l fc=2 . 



ci;mod\2 

I . 1 , / i x2 



(A3) 



In order to estimate the parameters that fit the data given the proposed model we have to maximise the likelihood or, 
equivalently, minimise ln£. Let p = [B\\, B±, 4>] denote the set of parameters of our model. Following the standard approach, 
to estimate the vector of parameters p, we must solve the following set of equations: 



a In £. _ g _ o 

"I ~ ~ ^ un 

1=1 j — 1 k—2 



dpi 



ni;mod\2 

- + ln(cr ifc ) 



V/. 



(A4) 



There is an important point to clarify given that, in our cas e, the model is not fully analytical because it depends on the 
observed intensity profile (|Asensio Ramos fc Manso Sainj|201ll ). Therefore, assuming no correlation between the model and 
the observable (the only difference being produced by uncorrelated Gaussian noise), the variance for each Stokes parameter 
is given by: 



a 2 (5 2 ) + a 2 (5 2 mod ) 



a 2 {S 2 ) + C 4 B 4 ± cos 2 20a 2 (X") 
a 2 (S 3 ) + a 2 (Sr d ) = o- 2 (5 3 ) + C 4 B 4 ± sin 2 20a 2 (X") 



M 2 
M 2 

M 2 = a 2 (5 4 ) + a 2 (5 4 mod ) = a 2 ^) + C 2 B\g 2 {X') 



(A5) 
(A6) 
(A7) 



Thanks to the previous equations, the variances that appear in Eq. (A4) depend on the actual parameters, which makes the 
minimisation much harder. Luckily, in the weak field regime and for the observational spectral resolution of interest nowadays, 
it is easy to verify that a 2 (5 l mod ) < a 2 (S z ). In any case, this condition should be checked before carrying out any inversion 
using the formulae derived in this work. Assuming a first order approximation to the derivative T' , we find: 



where K = 1 for the resolved case and K = 6-2u-3v ^ or ^ ne dipole case. Assuming a (S" 
in intensity is the same as in the circular polarisation, it holds that: 



Ax > V2CB }l KAg, 



(A8) 

^ a 2 (Si) and that the noise 



(A9) 



which means that the spectral sampling has to be larger than or of the order of the Zeeman splitting. For instance, for a 
wavelength of 5000 A, a Lande factor of 1.5, and Bu = 500 G, the spectral sampling has to be larger or equal to 12 mA 
for K = 1. Luckily, this is the case in most of the observational cases. For example, for the same wavelength, the expected 
sampling for two spectrographs with a resolving power of R — 60000, 300000 are 83 mA and 16 mA, respectively. For the 
linear polarisation we have, if we assume that their associated noise equals the noise in intensity, that: 

Ax 5s x/2GK' cos 2(j)CB±A 



where K' — 1/4 for the resolved case and K' = 
we end up with Ax = 6 mA. 



Ax ^ x/2GK' sin2(j>CB±A, (A10) 
i 420- 68^-105^ for the dipole case For B± = 500 G] ^ = ° and G = 2.25, 
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Taking the previous considerations into account, Eq. (A4) can be simplified to: 



dlnC 
dpi 



n t n x 4 



i=l j=l fc=2 



dpi 



{sL-s 



i;mod\ 2 
jk > 



HJ 2 



dpi ' 



(All) 



where the well known \ 2 merit function is defined as: 



x ^ K ? ) 2 z - 



/ 7 yi 7 yi;mod\2 



(A12) 



Explicitly, the derivatives of the \ 2 with respect to the parameters we want to infer are: 



i \2 



|£ = AC 2 B ± (cos 20 X; + sin 20 £ + AC^l ( cos 2 20 £ ^ + sin 2 20 £ 

5£ =4C 2 *?i f COS 20X^S - sin2 ^ES^) +4CXsin20cos20 ^ ^ 



(A13) 
(A14) 
(A15) 



By forcing these derivatives to zero, we obtain the maximum-likelihood estimate of each parameter. For the longitudinal 
magnetic field, using Eq. A13, we obtain a unique solution: 



(A16) 



Dividing Eq. A14 and A15, we obtain a solution for the azimuth: 



tan 20 



{i")) 2 U]Z") 
^ (-b,) 2 ^ < CT »,) 2 



(A17) 



Dividing Eq. A15 by cos 20 and using Eq. A17, after some algebra we obtain two solutions. One is B± = 0, which is not valid 
since this solution maximises the \ 2 - We have tested this computing the second derivative. The solution that minimises the 
X 2 is given by: 



^ < CT «,) 2 ^ ( CT b,) 2 / V 2 ^' ( CT b,) 2 ^ < CT «,) 2 



C 2 



^ ("hi) 2 ^ (" 



(A18) 



The errors associated with each parameter are computed assuming that the surface around the maximum likelihood is 
approximately a multidimensional Gaussian. This is equivalent to assuming a parabolic approximation to the \ 2 close to the 
minimum. The curvature close to the minimum is given by the Hessian matrix £, whose elements are: 



i a 2 x 2 

(,kl — 77- 



2 dpkdpi 



(A19) 
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The matrix can be calculated using the following second derivatives: 

o2 2 

d X 



8C 4 23 2 



0fijj 
dB\ 

d 2 x 2 

<90 2 

dx 2 

dB\\dB±_ 

dx 2 

dB\\dcf> 

dx 2 

dBxd<j> d<j>dBA 



2C 2 J2^^- 



8C 4 B 4 ± 

dx 2 

dB±dB l 

dx 2 

d<j>dB\\ 

dx 2 



sin 

= 
= 



2^E 



(4 



(1" 



(4 



+ cos' 



ij 

! 2« 



(I" 



(I" 



= 8C 4 Sl sin 20 cos 20 



(4 



(A20) 
(A21) 
(A22) 

(A23) 
(A24) 

(A25) 



The square root of the diagonal of the covariance matrix gives the error estimates for each parameter. Such covariance 
matrix is just the inverse of the Hessian matrix, C = C _1 - ln our case: 



C = 










C22 C23 
C32 C33 



with 



Cll = 



C22 = 



c 2 £ 



sin 2 20 E 



+ cos 2 20 £ 



(i"p 2 



4C 4 23 2 



cos 2 20 sin 2 20 E 



(i"j) 2 



E 



(i"P 2 
« (^) 2 



+ E 



sin 20 cos 20 I E 



C23 = C32 = 



C33 



-E 



(i"j) 2 
« (-b,) 2 



E 



2 



cos 2 20 sin 2 20 E 



: 20E 



E 



+ E 



y (-0,) 2 



E 



(I"') 2 



2 + sin 2 20 E 



(i"j) 2 



4C 4 23 4 



cos 2 20 sin 2 20 E 17 



E 



E 



(i"j) 2 



The error bars are thus given by: 



= ±AVC~ k 



(A26) 



(A27) 



(A28) 



(A29) 



(A30) 



(A31) 



with A = 1, 1.65, 2, 2.57, 3, and 3.89 for a confidence level of 68.3%, 90%, 95.4%, 99%, 99.73%, and 99.99%, respectively. This 
computation assumes that, to estimate the error of a parameter, we fix the values of the rest to the ones that m aximize the 
likeli hood and compute the confidence levels for the one-dimensional probability distribution of this parameter (see lPress et al.l 
1 19861 . for more details). Note that this approximation does not take into account the degeneracies between p arameters since 
it do es not integrate the probability distribution of the rest of the parameters but fixes a certain value (see lAsensio Ramosl 
120111 . for a robust Bayesian inversion). Note also that the covariance matrix is diagonal if the standard deviation of Q and 
1A are the same = &uj)- This is generally the case in solar observations in which the measurement efficiencies for linear 
polarization are similar. 



